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The  interaction  of  a  premixed  methane/ air  flame  with  a  counter-rotating  vortex  pair  is  analyzed  using  a 
parallel  low-Mach-number  computational  model  that  is  based  on  a  detailed  C4C2  chemical  mechanism. 
Attention  is  focused  on  the  transient  response  of  the  heat  release  rate  and  the  flame  structure  at  the 
centerline  of  the  vortex  pair.  Results  are  obtained  for  vortex  pairs  of  different  strengths  under  lean,  stoi¬ 
chiometric,  and  rich  conditions.  For  the  range  of  vortex  strengths  considered,  the  computations  indicate 
that  the  heat  release  rate  in  the  rich  flame  decays  significantly  faster  than  in  the  stoichiometric  flame;  this 
behavior  is  consistent  with  recent  experimental  measurements.  Meanwhile,  the  heat  release  rate  in  the 
lean  flame  decays  at  a  slightly  slower  rate  than  in  the  stoichiometric  flame.  The  transient  response  of  flame 
radicals  such  as  H,  CH,  OH,  and  HCO  is  also  analyzed.  The  analysis  reveals  a  complex  nonlinear  depen¬ 
dence  of  the  transient  structure  on  both  the  vortex  strength  and  the  stoichiometry. 


Introduction 

Unsteady  interactions  of  premixed  flames  with 
vortical  structures  are  an  essential  feature  of  many 
reacting  flow  applications,  including  burners,  com¬ 
bustors,  furnaces,  and  fires.  Consequently,  numer¬ 
ous  efforts  have  focused  on  isolating  these  interac¬ 
tions  and  on  analyzing  flame  response  to  unsteady 
strain. 

One  common  approach  to  addressing  this  issue  is 
based  on  one-dimensional  studies  of  opposed  jet 
flames  with  imposed  sinusoidal  perturbations  [1-3], 
These  studies  have  established  a  number  of  inter¬ 
esting  trends,  in  particular  concerning  the  effects  of 
imposed  strain  amplitude  and  frequency  and  of  the 
mixture  stoichiometry.  Dynamic  features  have  also 
been  identified,  including  extinction,  partial  extinc¬ 
tion,  and  reignition. 

A  second  approach  is  based  on  analyzing  flame 
interactions  with  vortices  [4-6],  counter-rotating 
vortex  pairs  [7—11],  or  toroidal  vortex  rings  [12-14], 
These  investigations  have  also  led  to  several  obser¬ 
vations  regarding  unsteady  heat  release  and  have  in 
addition  enabled  the  study  of  multidimensional  phe¬ 
nomena,  including  flame  stretch  and  contortion, 
pocket  formation,  and  baroclinic  vortieity  genera¬ 
tion. 

There  appears  to  be  a  strong  similarity  between 
the  one-dimensional  strained  flame  and  the  local 
flame  zone  at  the  centerline  of  a  counter-rotating 
vortex  pair.  In  both  situations,  the  flame  is  subjected 


to  an  unsteady  strain  field,  while  curvature  effects 
are  naturally  minimized.  Consequently,  one  may  be 
inclined  to  expect  a  similar  flame  response,  at  least 
in  a  qualitative  fashion.  However,  this  similarity  may 
not  necessarily  hold.  A  specific  example  was  pro¬ 
vided  in  our  recent  effort  [15],  which  focused  on  the 
interaction  of  rich  and  stoichiometric  flames  with  a 
counter-rotating  vortex  pair.  Consistent  with  exper¬ 
imented  results,  the  computations  in  Ref.  [15]  found 
that  the  rich  flame  exhibits  a  faster  response  to  the 
unsteady  strain  than  the  stoichiometric  flame.  This 
behavior  is  in  contrast  to  the  opposed-jet  flame  study 
in  Ref.  [2] ,  which  noted  that  the  sinusoidally  excited 
rich  flame  has  a  slower  response  than  the  stoichio¬ 
metric  flame. 

The  above  results  raise  a  variety  of  questions  re¬ 
garding  the  origin  of  the  observed  differences  be¬ 
tween  the  strained  one-dimensional  flame  results 
and  the  strained  flame  at  the  centerline  of  the  vortex 
pair.  One  possible  cause  behind  these  differences 
concerns  the  imposed  strain,  whose  large  values  in 
the  two-dimensional  flow  may  prevent  us  from  using 
a  linearized  analysis  to  establish  a  relationship  with 
the  one-dimensional  flame  under  sinusoidal  pertur¬ 
bations.  Another  possibility  concerns  the  flame  con¬ 
figuration  itself,  which  is  typically  anchored  in  the 
one-dimensional  setting  and  freely  propagating  in 
the  two-dimensional  case. 

Motivated  by  the  desire  to  resolve  and  reconcile 
these  differences,  the  present  work  focuses  on  the 
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effects  of  strain  amplitude  on  the  interaction  be¬ 
tween  a  premixed  methane/air  flame  and  a  counter¬ 
rotating  vortex  pair.  Parallel  computations  are  per¬ 
formed  using  a  two-dimensional,  low-Mach-number 
reacting  flow  model  that  is  based  on  a  detailed  C  |  - 
C2  chemical  mechanism.  As  further  discussed  below, 
we  consider  vortex  pairs  of  different  strengths,  and 
in  each  case  results  are  obtained  at  stoichiometric 
(0  =  1),  rich  (0  —  1.2),  and  lean  (0  =  0.8)  con¬ 
ditions.  The  flame  response  is  analyzed  by  contrast¬ 
ing  the  transient  evolution  of  the  heat  release  rate 
as  well  as  flame  radicals,  including  CH,  OH,  H,  and 
HCO.  These  detailed  computations  are  made  pos¬ 
sible  using  recent  algorithmic  developments,  which 
are  highlighted  in  the  following  section. 


Formulation  and  Numerical  Schemes 

The  computational  model  used  in  the  simulations 
below  is  adapted  from  our  previous  efforts  in  Refs. 
[16,17],  which  provide  a  detailed  discussion  of  the 
numerical  construction.  Brief  descriptions  of  the 
physical  formulation  and  numerical  scheme  are  pro¬ 
vided  below. 


Physical  Model 


The  physical  model  is  based  on  the  zero-Mach- 
number,  compressible  flow  approach  [18]  which  ig¬ 
nores  acoustic  waves  and  allows  us  to  decompose  the 
pressure  into  a  spatially  uniform  component  P0(f) 
and  a  hydrodynamic  component  p(x,t)  which  varies 
both  in  space  and  time.  We  assume  a  two-dimen¬ 
sional  open  domain,  a  gas  mixture  with  zero  bulk 
viscosity,  and  a  detailed  chemical  kinetic  mechanism 
involving  N  species  and  K  elementary  reactions.  So- 
ret  and  Dufour  effects  are  ignored,  as  are  body 
forces  and  radiant  heat  transfer.  The  mixture  is  also 
assumed  to  obey  the  perfect  gas  law,  with  individual 
species  molecular  weights,  specific  heats,  and  en¬ 
thalpies  of  formation. 

Under  the  above  assumptions,  the  governing 
equations  are  expressed  in  non-dimensional  form  as 
[16,17] 
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respectively.  Here,  p  is  the  density;  T  is  the  tem¬ 
perature;  v  =  (u,  v)  is  the  velocity  vector;  Y;  is  the 
mass  fraction  of  species  i;  p  is  the  dynamic  viscosity; 
X  is  the  thermal  conductivity;  cp  is  the  mixture  spe¬ 
cific  heat;  gj  is  the  diffusive  mass  flux  of  species  i ; 

N 

G  =  2  cp,ig; 


cp ;  is  the  specific  heat  of  species  i;  Wj  and  Wj  are  the 
molecular  weight  and  chemical  production  rate  of 
species  i,  respectively; 


w  =  i/  2  w* 


wT  is  the  rate  of  chemical  heat  release;  Re,  Pr,  Sc, 
and  Da  are  the  Reynolds,  Prandtl,  Schmidt,  and 
Damkohler  numbers,  respectively;  and  0X  and  0,, 
are  the  viscous  stress  terms.  For  computational  ef¬ 
ficiency,  and  given  that  the  flow  cases  under  study 
involve  significant  N2  dilution,  mixture  transport 
properties  are  set  to  those  of  N2,  and  species  diffu¬ 
sion  coefficients  are  set  to  those  into  N2  at  the  local 
temperature. 


Numerical  Scheme 

Numerical  simulation  of  the  governing  equations 
is  performed  using  a  stiff,  operator-split,  finite-dif¬ 
ference  scheme  whose  construction  is  detailed  in 
Ref.  [17],  The  scheme  is  a  variable-density  projec¬ 
tion  method  for  reacting  flow,  which  relies  on  a  sec¬ 
ond-order  spatial  discretization  of  the  governing 
equations.  Field  variables  are  discretized  using  a 
staggered  grid  with  uniform  cell  size  along  each  co¬ 
ordinate  direction.  Velocity  components  are  speci¬ 
fied  at  cell  edges,  while  scalar  variables  are  specified 
at  cell  centers. 

One  of  the  key  features  of  the  computational  con¬ 
struction  in  Ref.  [17]  is  the  use  of  a  symmetric  split¬ 
ting  of  the  diffusion  and  reaction  terms  and  the  in¬ 
corporation  of  a  stiff  scheme  for  the  numerical 
integration  of  the  reaction  source  term.  The  splitting 
is  implemented  in  such  a  way  that  the  diffusion 
terms  are  updated  in  two  half-time-step  integrations, 
which  are  separated  by  a  full-time-step  integration 
of  the  reaction  terms.  Each  diffusion  half-step  com¬ 
prises  several  fractional  substeps,  which  allows  the 
use  of  a  globed  time  step  that  is  several  times  larger 
than  the  critical  diffusional  time  step.  Meanwhile, 
the  reaction  source  terms  in  the  species  and  density 
equations  are  handled  using  a  stiff  integrator,  which 
is  adapted  from  the  DVODE  package  [19], 
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The  convergence  properties  and  computational 
performance  of  the  stiff,  split  scheme  are  also  dis¬ 
cussed  in  Ref.  [17].  In  particular,  tests  demonstrate 
that  the  construction  is  effective  in  defeating  the 
stiffness  of  both  the  diffusion  and  reaction  source 
terms.  This  leads  to  order-of-magnitude  speedup 
over  explicit,  non-stiff  schemes  and  enables  large- 
scale  unsteady  computations  of  reacting  flow  with 
detailed  kinetics. 

Results  and  Discussion 

Computations  of  the  interaction  of  a  premixed 
methane/air  flame  with  a  two-dimensional  vortex 
pair  are  performed  in  a  rectangular,  0.4  X  1.6  cm 
domain.  The  domain  is  discretized  using  a  uniform 
grid  with  256  X  1024  cells  and  mesh  size  h  — 
15.625  /an.  This  fine  resolution  level  is  needed  to 
adequately  resolve  the  flame  structure.  The  initial 
condition  consists  of  the  superposition  of  the  velocity 
field  induced  by  a  periodic  row  of  vortex  pairs,  with 
temperature,  density,  and  mass  fraction  distributions 
of  a  steady-propagating,  one-dimensional  premixed 
flame  computed  with  CHEMKIN  [20,21],  We  use 
the  GRIMech  1.2  [22]  mechanism  with  32  species 
and  177  reactions.  The  flame  is  initially  flat  and  hor¬ 
izontal,  with  reactants  below  and  products  on  top. 
The  vortex  pair,  initially  in  the  reactants,  propagates 
upward  into  the  flame. 

We  consider  N2-diluted  methane/air  reactants 
mixtures  with  varying  stoichiometry.  Included  in  the 
analysis  are  lean  (0  =  0.8),  rich  (0  =  1.2),  and 
stoichiometric  mixtures.  In  all  cases,  a  20%  N2  di¬ 
lution  level  is  used.  In  addition  to  the  stoichiometry, 
we  vaiy  the  initial  characteristics  of  the  vortex  pairs. 
In  all  cases,  the  initial  distances  between  the  vortex 
centers  in  a  pair  is  held  fixed,  3  =  0.25  cm.  On  the 
other  hand,  the  strengths  of  the  vortices  are  varied 
by  up  to  a  factor  of  10;  we  use  r  =  0.07,  0.028,  and 
0.007  nr/s.  Since  the  self-induced  velocity  of  the 
vortex  pair  is  proportional  to  its  strength,  these  cases 
are  referred  to  as  fast,  moderate,  and  slow,  and  la¬ 
beled  F,  M,  and  S,  respectively.  As  will  be  shown 
later,  the  peak  strain  rate  field  experienced  by  the 
flame  also  scales  proportionally  with  the  vortex 
strength. 

Flame-Vortex  Interaction 

In  all  cases,  simulations  are  carried  over  a  0.5  ms 
period  in  order  to  observe  the  transient  flame  re¬ 
sponse  to  the  strain  rate  field  induced  by  the  vortex 
pairs.  Fig.  1  shows  the  heat  release  rate  and  vorticity 
fields  at  the  end  of  the  computations  for  all  nine 
cases.  The  interaction  between  the  flame  and  the 
vortex  pair  leads  to  the  generation  of  opposite  sign 
vorticity  [4,10],  which  is  depicted  in  Fig.  1  using 
dashed  contours.  Note  that  identical  contour  (vor¬ 
ticity)  and  shading  (temperature)  levels  are  selected 


in  all  frames.  As  the  strength  of  the  vortex  pair  is 
increased,  the  number  of  vorticity  contours  in¬ 
creases;  the  vortex  penetration  distance  into  the 
flame  and  the  extent  of  flame  stretching  increase  as 
well.  Note  that  in  case  S,  the  interaction  is  still  in  its 
early  stages,  flame  deformation  is  weak,  and  a  robust 
flame  region  is  observed  at  all  values  of  0.  For  case 
F,  however,  significant  contortion  of  the  flame  oc¬ 
curs,  and  it  is  evident  from  Figure  1  that  the  peak 
heat  release  rate  in  the  rich  flame  has  dropped  sig¬ 
nificantly.  This  indicates  that  flame  response  de¬ 
pends  substantially  on  both  the  vortex  strength  and 
the  stoichiometry. 

Unsteady  Strain  Rate 

In  order  to  analyze  these  changes,  we  first  quantify 
the  unsteady  strain  rate  that  the  flame  experiences 
during  the  interaction.  We  focus  our  attention  on  the 
centerline  of  the  vortex  pair,  where  the  stretching  of 
the  flame  is  large  and  the  curvature  of  the  front  is 
small.  We  choose  the  location  of  10%  Yfuei_max  as  a 
reference  flame  location,  where  the  tangential  strain 
rate  is  computed.  The  results,  plotted  in  Fig.  2,  in¬ 
dicate  that  the  flame  is  subjected  to  strain  rates  that 
increase  in  proportion  to  increased  vortex  circulation 
in  cases  S,  M,  and  F.  For  a  given  initial  vortex 
strength,  the  strain  rates  vary  with  time  and  stoichi¬ 
ometry.  This  is  most  accentuated  in  the  F  cases, 
where  significant  flame  contortion  and  baroclinic 
vorticity  generation  occur  (see  Fig.  1).  In  these 
cases,  the  strain  rates  in  the  lean  and  rich  flames  are 
higher  than  in  the  stoichiometric  flame.  Results  sug¬ 
gest  that  this  development  is  a  consequence  of  the 
flame  interaction  with  the  flow.  Earoclinic  vorticity 
generation  at  the  flame  leads  to  the  formation  of  a 
vortex  dipole  in  the  vicinity  of  the  primary  vortex 
[4,10],  The  results  suggest  that  the  strength  of  the 
opposite-sign  vortex  dipole  and  its  position  during 
the  interaction  result  in  the  observed  difference  in 
the  strain  rate  levels. 

On  the  other  hand,  the  variations  of  strain  with 
time  and  stoichiometry  are  relatively  weak  and  fall 
within  ±  10%  of  the  mean.  This  provides  a  reason¬ 
able  basis  for  examining  the  effect  of  stoichiometry 
on  the  transient  flame  response,  as  the  strain  rate  is 
roughly  the  same  for  the  three  stoichiometries  in 
each  of  the  F,  M,  and  S  cases. 

Heat  Release 

The  global  flame  response  is  interpreted  in  terms 
of  the  peak  heat  release  rate,  whose  evolution  is  plot¬ 
ted  in  Fig.  3.  The  profiles  are  normalized  by  their 
initial  values  in  each  case.  The  results  indicate  that 
the  peak  heat  release  rate  decays  as  the  vortex  pair 
impinges  into  the  flame.  Consistent  with  the  trend 
in  the  strain  rates,  the  rate  of  decay  in  heat  release 
is  faster  as  the  vortex  strength  increases.  Fig.  3  also 
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Fig.  1.  Heat  release  rate  (shading)  and  vorticity  (con¬ 
tours)  distributions  at  t  =  0.5  ms.  Shown  are  results  at  lean 
(left),  stoichiometric  (center),  and  rich  (right)  conditions, 
for  strong  (top),  intermediate  (middle),  and  weak  (bottom) 
vortices. 


Time  (ms) 


Fig.  2.  Tangential  strain  rate  time  history  at  the  flame 
10%  Tfuei>max  location  on  the  vortex  pair  centerline.  Curves 
correspond  to  the  slow  vortex  pair  at  lean  (S0.8),  stoichio¬ 
metric  (S1.0)  and  rich  conditions  (SI. 2);  the  intermediate 
vortex  pair  at  lean  (MO. 8),  stoichiometric  (M1.0)  and  rich 
conditions  (M1.2);  and  the  fast  vortex  pair  at  lean  (F0.8), 
stoichiometric  (F1.0),  and  rich  conditions  (FI. 2). 


Fig.  3.  Evolution  of  the  normalized  peak  heat  release 
rate  on  the  vortex  pair  centerline.  Curves  are  labeled  as 
in  Fig.  2. 


shows  that,  in  each  flow  case,  the  rich  flame  responds 
significantly  faster  than  either  the  lean  or  stoichio¬ 
metric  flames.  In  earlier  work  [15]  we  noted  that  for 
case  F  the  heat  release  in  the  rich  flame  responds 
faster  than  in  the  stoichiometric  flame,  in  contrast  to 
earlier  one-dimensional  flame  results  [2],  The  pres¬ 
ent  computations  enable  us  to  amplify  the  previous 
finding,  since  the  same  observation  holds  even  as  the 
strain  rate  is  varied  (reduced)  over  a  tenfold  range. 
Remarkably,  the  present  computations  also  indicate 


that  the  lean  flame,  which  is  less  robust  than  the 
stoichiometric  flame,  exhibits  slower  heat  release  de¬ 
cay  rate  than  the  latter.  This  observation  also  holds 
for  all  the  flow  speeds  considered. 

We  have  also  suggested  in  Ref.  [15]  that  the  global 
flame  response  time  is  governed  by  the  relative  roles 
of  H  and  OH  in  the  kinetics  and  the  faster  response 
time  of  H  versus  OH.  The  lean  and  rich  flames  are 
both  weak  relative  to  the  stoichiometric  flame.  At  t 
=  0,  peak  heat  release  rates  are  147,  306,  and  176 
cal/cm3-s  for  </>  =  0.8,  1.0,  and  1.2,  respectively. 
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Time  (ms) 

Fig.  4.  Evolution  of  the  normalized  peak  H  mole  frac¬ 
tion  on  the  centerline  of  the  vortex  pair.  Curves  are  labeled 
as  in  Fig.  2. 


Fig.  5.  Evolution  of  the  normalized  peak  OH  mole  frac¬ 
tion  on  the  centerline  of  the  vortex  pair.  Curves  are  labeled 
as  in  Fig.  2. 


Thus,  robustness  (e.g.,  peak  fuel  consumption  or 
heat  release  rates)  would  not  explain  the  trends  in 
the  figure.  On  the  other  hand,  we  find  that  for  hy¬ 
drocarbon  species  on  the  path  to  HCO  (thus  asso¬ 
ciated  with  primary  heat  release  in  the  fuel  con¬ 
sumption  layer),  for  example,  CH4,  CH20,  and 
HCO,  the  ratio  of  carbon  flux  carried  by  H  reactions 
to  that  carried  by  OH  reactions  increases  with  <t>. 
This  suggests  that  the  relative  role  of  H  versus  OH 
in  the  fuel-consumption  chemistry  is  significantly 
enhanced  in  going  from  lean  to  stoichiometric  to 
rich,  consistent  with  the  change  in  response  time  of 
the  heat  release.  Moreover,  the  ratio  of  available  H 
to  OH  increases  accordingly.  Thus,  since  H  has  a 


faster  response  than  OH  in  this  flame,  the  heat  re¬ 
lease  rate  responds  at  faster  timescales  as  <P  in¬ 
creases,  due  to  the  decreased  role  of  OH  and  in¬ 
creased  role  of  H.  Also  note  that  the  pathway  from 
CH3  to  CH20  relies  primarily  on  O,  rather  than  H 
and  OH.  On  the  other  hand,  we  find  the  response 
of  O  to  be  as  fast  or  faster  than  that  of  H  in  all  cases, 
and  so  O  does  not  appear  to  play  a  limiting  role  in 
the  response  of  the  heat  release  rate. 

OH  and  H 

The  response  of  H  and  OH,  alluded  to  above,  is 
illustrated  in  Figs.  4  and  5.  The  results  show  that  the 
OH  response  rate  approaches  that  of  II  only  at  0  = 
1.2,  but  is  still  slower  even  then.  The  results  also 
indicate  that  the  response  of  H  and  OII  have  com¬ 
plex  dependence  on  d>  and  r,  which  does  not  coin¬ 
cide  with  the  response  of  the  heat  release  rate.  Spe¬ 
cifically,  while  H  mole  fraction  decays  fastest  at  (I> 
=  1.2  in  cases  S  and  M,  in  case  F  it  is  the  lean  flame 
that  exhibits  fastest  H  decay.  This  indicates  that  in¬ 
teresting  changes  occur  at  large  strain  amplitudes. 

It  is  also  interesting  to  note  that  the  stoichiometric 
flame  has  the  slowest  H  decay  rate,  in  contrast  to 
the  above  observations  for  wT.  Meanwhile,  Fig.  5 
shows  that  at  fixed  /  ,  OH  decays  fastest  in  the  rich 
flame,  in  agreement  with  the  heat  release;  this  sup¬ 
ports  the  observation  that  OH  is  a  limiting  factor  in 
overall  flame  response.  The  lean  and  stoichiometric 
flames,  on  the  other  hand,  have  similar  OH  response 
in  cases  S  and  M,  while  the  lean  case  is  slightly  faster 
in  case  F ;  this  behavior  is  not  consistent  with  that  of 
the  heat  release  rate. 

CHoO,  HCO ,  and  [OH][CHsO] 

The  response  of  peak  CHzO  (not  shown)  is  similar 
to  that  of  OH  and  exhibits  similar  dependence  on  </> 
and  r.  On  the  other  hand,  the  response  of  peak 
HCO  (also  omitted)  closely  coincides  with  that  of 
heat  release  rate  (Fig.  3).  This  is  in  agreement  with 
and  extends  the  correlation  established  between 
HCO  mole  fraction  and  heat  release  rate  in  Refs. 
[11,23],  As  outlined  in  Ref.  [9],  experimental  chal¬ 
lenges  associated  with  measurement  of  HCO  moti¬ 
vate  the  use  of  the  peak  concentration  product  of 
[OH][CHaO]  as  surrogate.  This  product  was  in  fact 
found  to  be  an  acceptable  surrogate  of  both  HCO 
and  icx  in  [9,23] .  The  same  trend  is  observed  in  the 
present  computations  at  all  values  of  d>  and  r,  which 
amplifies  the  scope  of  previous  findings. 

ch3o 

The  complexity  of  the  flame  response  is  also  re¬ 
flected  in  the  evolution  of  the  methoxy  radical, 
CH30,  shown  in  Fig.  6.  CH3O  exhibits  transient  ac¬ 
cumulation  whose  amplitude  increases  with  71  For 
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Fig.  6.  Evolution  of  the  normalized  peak  CH30  mole 
fraction  on  the  centerline  of  the  vortex  pair.  Curves  are 
labeled  as  in  Fig.  2. 


Time  (ms) 


Fig.  7.  Evolution  of  the  normalized  peak  CO  mole  frac¬ 
tion  on  the  centerline  of  the  vortex  pair.  Curves  are  labeled 
as  in  Fig.  2. 

the  lean  F  case,  a  fivefold  increase  in  peak  CH30 
mole  fraction  is  achieved  even  as  the  flame  proceeds 
toward  extinction.  At  fixed  r,  the  methoxy  mole  frac¬ 
tion  exhibits  shorter  rise  times  in  the  lean  and  rich 
cases  than  in  the  stoichiometric  flame.  On  the  other 
hand,  it  achieves  highest  amplitude  of  transient  ac¬ 
cumulation  in  the  lean  case,  which  has  the  smallest 
peak  heat  release  rate.  The  more  robust  stoichio¬ 
metric  flame  is  the  one  having  the  smallest  ampli¬ 
tudes  and  the  longest  rise  times. 

CH 

The  peak  CH  mole  fraction  exhibits  a  similar  de¬ 
pendence  on  0  and  f  as  OH  (Fig.  5);  the  results  for 


CH  are  therefore  omitted.  In  particular,  as  for  icT, 
CH  decay  is  fastest  in  the  rich  flame.  The  response 
of  CH  in  the  lean  and  stoichiometric  cases  is  almost 
identical  in  cases  S  and  M.  On  the  other  hand,  as 
with  OH,  the  lean  case  exhibits  faster  decay  rate  of 
peak  CH  in  case  F;  this  behavior  is  not  in  agreement 
with  the  trends  established  for  wT. 

CO 

The  response  of  peak  CO  is  shown  in  Fig.  7,  which 
indicates  a  monotonous  decrease  in  speed  of  re¬ 
sponse  with  increasing  0.  This  behavior,  which  is 
observed  at  all  values  of  /  ',  is  in  contrast  with  the 
response  of  wT  (Fig.  3)  and  HCO.  This  is  somewhat 
surprising,  since  HCO  is  the  major  source  of  CO.  In 
fact,  we  do  find  that  the  response  of  the  peak  CO 
production  rate  increases  with  0  in  a  manner  similar 
to  HCO  and  wT.  This  suggests  that  the  cause  of  the 
inverse  behavior  of  CO  mole  fraction  is  due  to  trans¬ 
port,  rather  than  reaction.  Inspection  of  CO  mole 
fraction  profiles  reveals  that  the  amount  of  CO  be¬ 
hind  the  flame  increases  with  0,  dictated  by  the  cor¬ 
responding  increase  in  the  equilibrium  concentra¬ 
tion  of  CO  in  the  products.  The  absence  of  radiant 
heat  losses  in  the  model  leads  to  a  higher  plateau  of 
temperature  and  CO  in  the  products  region,  as  com¬ 
pared  with  a  non-adiabatic  flame  model.  Moreover, 
higher  levels  of  CO  are  obtained  at  higher  0.  This 
increase  in  the  available  CO  behind  the  flame  then 
leads  to  a  reduction  in  the  decay  rate  of  peak  CO 
mole  fraction  as  0  increases.  Even  with  radiant 
losses  included,  the  products  CO  level  in  the  im¬ 
mediate  vicinity  of  the  flame  is  expected  to  increase 
with  0 ,  so  the  role  of  transport  in  the  decay  rate  of 
CO  may  also  be  observed  in  a  non-adiabatic  setting. 
This  is  a  subject  for  future  work. 

Conclusions 

A  parametric  study  is  conducted  of  the  interaction 
of  a  premixed  methane/ air  flame  with  a  counter-ro¬ 
tating  vortex  pair  under  a  range  of  reactants  stoichi¬ 
ometry  and  flow  time  scales.  Attention  is  focused  on 
the  transient  response  of  the  heat  release  and  flame 
structure  on  the  vortex  pair  centerline.  In  particular, 
computed  results  indicate  the  following. 

1.  For  all  values  of  r  considered,  the  heat  release 
timescale  decreases  as  0  increases.  This  trend 
generalizes  earlier  findings  in  Ref.  [15],  which 
were  restricted  to  fast  vortices  and  to  stoichio¬ 
metric  and  rich  conditions.  Remarkably,  and  in 
contrast  with  results  for  one-dimensional  flames 
under  sinusoidal  stretching,  the  slower  propagat¬ 
ing  rich  flame  responds  faster  than  the  stoichio¬ 
metric  flame.  Meanwhile,  the  lean  flame  re¬ 
sponds  at  a  slightly  slower  rate  than  the 
stoichiometric  flame.  Consequently,  the  response 
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time  scale  does  not  generally  correlate  with  peak 
heat  release  rate  or  the  fuel  consumption  rate. 

2.  The  response  of  the  heat  release  rate  appears  to 
be  governed  by  the  relative  roles  of  H  and  OH 
in  the  kinetics.  As  0  increases,  more  of  the  fuel 
consumption  pathways  are  dominated  by  H  re¬ 
actions  due  to  the  scarcity  of  OH  and  therefore 
respond  at  the  faster  H  time  scales. 

3.  At  all  values  of  0  and  F,  the  peak  HCO  mole 
fraction  exhibits  a  response  similar  to  that  of  the 
peak  heat  release  rate.  By  including  a  wide  range 
of  flow  timescales  as  well  as  lean  conditions,  the 
present  results  extend  the  scope  of  earlier  cor¬ 
relations  [11,23], 

4.  The  responses  of  H  and  OH  exhibit  some  simi¬ 
larities  with  that  of  heat  release  rate  but  reveal 
differences  which  depend  on  71  Meanwhile,  the 
response  of  CH30  reveals  transient  accumula¬ 
tions  with  highest  amplitude  occurring  in  the  lean 
flame.  The  decay  of  peak  CO  is  fastest  in  the  lean 
flame,  in  contrast  with  wT,  HCO,  and  the  CO 
production  rate.  The  results  suggest  that  trans¬ 
port  plays  a  dominant  role  in  the  observed  CO 
behavior. 

The  observed  dependence  of  the  flame  response 
on  strain  and  stoichiometry  does  not  appear  to  follow 
a  simple  correlation  based  on  global  parameters, 
such  as  heat  release  rate,  propagation  speed,  or 
flame  thickness.  In  particular,  the  computations  in¬ 
dicate  that  the  flame  kinetics  play  a  dominant  role 
in  the  transient  behavior  of  the  heat  release  rate. 
This  prevents  us  from  adapting  simplified  interpre¬ 
tations  based  on  diffusion  timescales.  Despite  con¬ 
sidering  a  wide  range  of  flow  timescales  and  varying 
the  mixture  stoichiometry,  a  sharp  contrast  exists  be¬ 
tween  the  present  two-dimensional  computations 
and  earlier  results  for  one-dimensional  flames  under 
sinusoidal  forcing.  The  present  experiences  suggest 
that  this  discrepancy  is  unlikely  to  be  due  to  the 
strain  rate  amplitudes.  An  extended  investigation  is 
currently  being  planned  to  further  explore  this  issue. 
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